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ABSTRACT 


The fuel temperature and reactor power responses following fuel -emitter separation 
were calculated using a one-dimensional, linearized, time -dependent -heat -transfer 
model. The time -dependent power density in the fuel is calculated using space - 
independent, monoenergetic neutron diffusion theory. Six neutron delay groups and 
temperature coefficients of reactivity for all diode components are included in the analy- 
sis. Fuel melting is predicted in the reference diode design, and radial cooling fins are 
recommended as a design modification to avoiding fuel melting. 
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SUMMARY 

Several thermionic reactor concepts have been investigated for electric -power gen- 
eration in space. One concept provides externally UOg -fueled thermionic diodes which 
are cooled by liquid lithium flowing in a coaxial channel at the center of the diode. Be- 
cause of the difference between the coefficients of thermal expansion for the fuel and 
emitter, temperature increases within the fuel might cause the fuel to "pull away" from 
the emitter. This situation might result in fuel melting and a subsequent nuclear catas- 
trophy. 

The consequences of fuel-emitter separation were analyzed mathematically by per- 
forming linearized transient heat transfer and neutronic calculations on a reference de- 
sign externally fueled thermionic diode. The purpose of the study is to determine if fuel 
melting would occur and how long it might take to initiate fuel melting. The results of 
the study indicate that the problem is severe enough to cause fuel melting. If fuel- 
emitter separation occurs in all the diodes in the reactor, the power rise accompany- 
ing separation could be detected by neutron flux monitors and reactor power could be re- 
duced in time to avoid fuel melting. If, however, fuel-emitter separation occurs in only 
a few diodes, reactor power changes will be too small to be detected by neutron flux 
monitors and melting of the fuel will occur until sufficient time elapses to allow the fuel 
to be vaporized and redeposited on the tungsten emitter. 

One method of avoiding fuel melting involves placing radial tungsten fins in the fuel 
region to allow conduction of heat from the fuel, through the fin, and into the emitter. 

If eight 6 -mil- (0. 015 cm) wide fins are incorporated in the reference diode design, fuel 
melting could be avoided after fuel-emitter separation. The redesign would also require 
3 percent increase in reactor volume and a 1.1 percent increase in fuel outer radius in 
order to maintain reactor criticality. 



INTRODUCTION 


Many power -source concepts are currently being investigated for electric power 
generation in space. Several of these concepts combine thermionic diodes which convert 
thermal energy to electrical energy and a nuclear power source which provides thermal 
energy from the fission of nuclear fuel. One thermionic reactor design discussed in ref- 
erence 1 incorporates a fueled -diode arrangement called an ” externally fueled flashlight” 
configuration. In this design, the fuel assemblies are hollow cylinders which surround a 
diode configuration with a central, coaxial cooling passage as shown in figure 1. The 
UOg fuel surrounding the diode is the source of thermal power. 

Ward and Ruch (ref. 2) recognized that the fuel will tend to separate from the diode 
when the temperature rises because the coefficient of thermal expansion of the fuel is 
higher than that of the diode emitter. If the fuel separates from the emitter, the resis- 
tance to heat conduction will increase and the fuel temperature will rise until the heat 
generated within the fuel can be transfered to the emitter. The temperature reached 
under these circumstances may be high enough to melt the fuel. Fuel melting must be 
avoided because it might result in damage to the clad due to rapid release of gaseous 
fission products; a greater loss of UO 2 fuel from fully vented fuel elements; a chemical 
reaction between the fuel and the clad causing a loss of clad strength; and a nuclear ex- 
cursion caused by redistribution of the fuel. 

The purpose of this study is to determine 

(1) The fuel temperature response as a result of fuel-emitter separation 

(2) How long it takes to initiate fuel melting (This information is required for the 

design of reactor safety circuits. ) 

(3) What design modifications might be made to the fuel assembly or in operating 

conditions either to prevent fuel-emitter separation or to prevent fuel melting. 

The reference design of the thermionic fuel assembly is described first. This is 
followed by a description of the mechanisms associated with fuel-emitter separation. 
Calculated values of the temperature response following separation are described next 
and various methods of avoiding fuel melting are discussed. 


SYMBOLS 


A area 

C reactivity coefficient 

AC thickness of fuel clad 

C p heat capacity 
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F mass flow rate 

f-p volume fraction of fuel in reactor 
h forced convection film coefficient 

k thermal conductivity 

L neutron lifetime 

l length of fuel cell 

M molecular weight 

m mass velocity of fuel vaporized 

n number of heat conducting fins 

P pressure 

Q/A heat flux 

q power density 

R specific radial field point 

AR spacing between fuel and emitter after accident 

r general radial field point 

T temperature 

t time 

U overall heat-transfer coefficient 

V volume 

w fin width 

$ total delayed neutron fraction 

e thermal emissivity 

T] diode efficiency 

r reactivity 

X decay constant 

p density 

a Stefan -Boltzmann constant 

6 angle 
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Subscripts: 

B coolant 

BE end clad -to -emitter 

BF end clad -to -fuel 

BP end clad 

b boundary between fuel and cooling fin 
C collector 

E emitter 

F fuel 

f initial value 

I insulator 

in coolant inlet 

max maximum 

o coolant outlet 

P Pipe 

r reactor core 

s surface 

DESCRIPTION OF REFERENCE DESIGN THERMIONIC FUEL ASSEMBLY 

A preliminary design of a fuel element containing externally fueled thermionic 
diodes was presented by Yacobucci (ref. 1) and a conceptual drawing of the fuel element 
is shown in figure 1. These fuel elements are uniformly distributed throughout a 
nuclear reactor and are electrically connected such that the desired electrical power 
output can be obtained. The diodes in each fuel element are stacked vertically along the 
length of the fuel element as illustrated in figure 1. The diodes are electrically con- 
nected in series and are physically separated by an electrical insulator. A typical diode 
contains a central coolant tube surrounded by the following concentric components: 

(1) An AlgOg electrical insulator which separates the collector from the coolant tube 

(2) A niobium collector 

(3) An interelectrode gap 

(4) A tungsten emitter 
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Figure 2. - Schematic drawing of externally 
fueled thermionic diode. 
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TABLE I. - NOMINAL CELL DIMENSIONS 
[Cell length = 1 inch (2. 54 cm).] 


Component 

Outer radius of 
component 

Material 

in. 

cm 

Coolant channel (Rj) 

0. 1605 

0.4077 

Li 

Pipe (R 2 ) 

.2005 

.5093 

Nb 

Insulator (R^) 

.2105 

. 5347 

A1 2°3 

Collector (R^) 

. 2360 

.5994 

Nb 

Interelectrode gap (R g ) 

.2460 

.6248 


Emitter (Rg) 

. 2660 

.6756 

w 

Fuel (R 7 ) 

.4290 

1.090 

uo 2 

Clad (R 8 ) 

.4490 

1. 140 

w 


TABLE H. - REACTOR OPERATING CONDITIONS 


Thermal 


Power density, Btu/(sec)(ft 2 ); W/cm^ 

3800; 142 

Coolant flow, lbm/sec; kg/sec 

0.079; O.O^ 

Emitter heat flux, Btu/(sec)(ft 2 ); W/cm 2 

68; 77 

Diode efficiency, percent 

11 

Nuclear 


Initial reactivity, F(0) 

0.0 

Neutron lifetime, sec 

4. OxlO -6 

% 

0. 22X10" 3 

* 2 

1. SOxlO”* 3 

q 

$3 

1. 33xl0” d 

q 

h 

2. 69xlO” J 

q 

h 

0. 78X10" J 


0. 28X10' 3 

^eff . 

Xp sec 1 

X 2 , sec" 1 
X^, sec 
X 4 , sec” 1 
X 5 , sec” 1 
Xg, sec” 1 

6.8X10' 3 
1. 24X10" 2 
3. 05X10 -2 
1. llxlO -1 
3.01X10 -1 
1. 13 
3.0 
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(5) Uranium dioxide fuel 

(6) Tungsten clad 

The components and dimensional symbols of a typical diode are schematically illus - 
trated in figure 2. The nominal diode dimensions used in the thermal analysis are pre- 
sented in table I. Reactor operating conditions are summarized in table n, along with 
some neutron kinetic properties of the reference reactor. 

DESCRIPTION OF FUEL-EMITTER SEPARATION MECHANISMS 

The mechanisms associated with the fuel-emitter separation problem may be defined 
as 

(1) The physical separation of the fuel from the emitter caused by some thermal 

perturbation in the fuel or emitter 

(2) Subsequent impairment of heat transfer from the fuel causing the fuel tempera- 

ture to increase and the emitter temperature to decrease 

(3) Possible reactor power excursion caused by the reactivity addition associated with 
a reduction in emitter temperature (negative temperature coefficient of reactivity) 
or an increase in fuel temperature (positive temperature coefficient of reactivity) 

Any perturbation in reactor power or coolant flow tending to increase the fuel tem- 
perature might cause the fuel to separate from the emitter because the coefficient of 
thermal expansion for the uranium dioxide fuel is about twice as large as that of the 
tungsten emitter. Therefore, the fuel expansion per degree temperature rise is greater 
than the emitter expansion and the fuel will tend to ’’pull away'* from the emitter when 
the fuel temperature increases. The bond between the fuel and emitter may not be 
strong enough to keep the fuel from separating from the emitter. After separation, the 
resistance to heat transfer between the fuel and emitter will increase. The size of the 
gap between the fuel and emitter is presented in figure 3 as a function of the rise in the 
fuel and emitter temperature. The gap size was calculated with coefficients of thermal 
expansion of 6. 74X10' 6 °R _1 (1. 22x10" 5 K" 1 ) and 3. 16X10" 6 0 R _1 (5. 68xl0" 6 K" 1 ) for 
the fuel (ref. 3) and emitter (ref. 4), respectively. 

The heat generated in the fuel can be dissipated by radiation to the emitter, mass 
transfer as a result of the fuel vaporizing and condensing on the cooler emitter, and 
conduction around the fuel-emitter gap through the end clad shown in figure 2. The fuel 
temperature will rise until the heat can be dissipated. The temperature rise may be 
sufficiently large to cause fuel melting. Fuel melting should be avoided because of the * 
potential danger from rapid release of fission products and chemical reaction between 
the fuel and clad. Also, if a sufficient number of fuel assemblies melt, the redistribu- 


7 




Fuel temperature rise, Tp - Tp 0 , °R 



Fuel temperature rise, Tp -Tp 0 , K 
Figure 3. - Gap formed when fuel separates from emitter. 

tion of fuel might introduce sufficient positive reactivity to cause a reactor power excur- 
sion and subsequent violent disassembly of the reactor structure. 

If the fuel emitter separation occurs in only a few assemblies, no detectable change 
in reactor power will occur. If, however, separation occurs in many fuel assemblies, 
reactor power will be affected because of the temperature coefficient of reactivity asso- 
ciated with the fuel assembly components. If the reactivity of the reactor increases as 
a result of the temperature change, reactor power and fuel temperature will increase. 
Because of the potential effect of reactivity coefficients on reactor power and tempera- 
ture, it is necessary to investigate the time -dependent neutronics response as well as 
the temperature response of the reactor. 

TEMPERATURE RESPONSE AFTER SEPARATION 

The temperature response of the fuel assembly components after fuel-emitter sepa- 
ration were calculated using the mathematical model described in appendix A. This 
model'allows the calculation of the one -dimensional (radial), time -dependent tempera- 
tures of each of the major components of the diode. Axial heat loss from the fuel to the 
end clad is permitted. All modes of heat transfer are approximated as linear functions 
of temperature and no mass transfer between the fuel and emitter is permitted. The 
diode emitter current density is assumed to be constant after separation. This is a 
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valid assumption because the diodes are electrically connected in series and, therefore, 
all experience the same current. The time -dependent power density in the fuel is cal- 
culated using space -independent, monoenergetic, neutron kinetics equations. Six neutron 
delay groups are included and the effect of component temperature variations on reactiv- 
ity (Doppler coefficients) are incorporated in the kinetics equation. 

The temperature and power responses were calculated for three separate conditions: 

(1) Fuel-emitter separation occurs in a single diode such that there is no appreciable 
change in reactor reactivity as the temperature of the diode components vary. 

(2) Fuel-emitter separation occurs in all diodes in the reactor and reactor reactiv- 
ity depends upon the temperature coefficients of reactivity and the temperature of the 
diode components. The temperature coefficients of reactivity for the reference design 
diode have not been calculated. 

The temperature coefficients of reactivity shown in the following table are compara- 
ble to those used by other investigators (unpublished information obtained by C. R. Fisher, 
A. J. Gietzen, C. A. Heath, W. G. Homeyer, and D. R. Mathews of General Atomics) 
for their study of the unit cell thermionic reactor design. 


Component 

Reactivity coefficient 

$/°R 

$/K 

Fuel 

Emitter 

Collector 

Insulator 

Pipe 

Coolant 

i 

0 

-1.6X10" 3 4 
-4. 1X10" 4 
-5. 3x10" 5 
-4. 1X10" 4 
0 

0 

-2. 9X10" 4 
-7. 4X10“ 4 
-9. 5xl0“ 5 
-7. 4xlO" 4 
0 


(3) Same as condition 2 except that the temperature coefficient of reactivity of the fuel 
is 8.0X10" 5 $/°R (1.4X10 -4 $/K). 

The initial (before accident) component temperatures are summarized in table HI. 
These temperatures were obtained from reference 1. The fuel -temperature rise, the 

emitter -temperature rise, and the relative reactor power are shown in figure 4 as a 
function of time after fuel-emitter separation. The reactor power is expressed relative 
to its initial value. The temperature changes in the coolant, coolant tube, insulator, 
and collector are much less than that of the fuel and emitter. Although these tempera- 
ture changes were included in the calculations, they are generally too small to be illus- 
trated in figure 4. 
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TABLE m. - INITIAL DIODE 
TEMPERATURES 


Temperature 

°R 

K 

Maximum fuel 

4860 

2700 

Emitter 

3600 

2000 

Collector 

2520 

1400 


Temperature coefficient of reactivity, $/°R 
Fuel Emitter Collector Insulator Pipe Coolant - 


0 0 0 0 0 0 

0 -1.6xl0" 4 -4.1xl0' 4 -5.3x10 -4. lxlD" 4 0 

8. 0x10 " 5 -1. 6xl0~ 4 -4,lxl0“ 4 -5.3XHT 5 -4. lxlO" 4 0 



Figure 4. - Temperature and power response when fuel separates from emitter in an externally fueled thermionic 
diode. 
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Following fuel-emitter separation, a thermal resistance is introduced between the 
fuel and emitter. Therefore, the heat transfer to all components located inside the fuel 
is reduced and the temperature of these components decreases. This is illustrated in 
all three cases in figure 4 by the 600° to 700° R (330 to 390 K) reduction in emitter tem- 
perature which occurs within 5 seconds after the accident is initiated. During the same 
period, the fuel temperature rises because the resistance to heat removal from the fuel 
is increased when fuel-emitter separation occurs. Because there is no provision in the 
mathematical model for inclusion of the latent heat of fusion, melting cannot be accounted 
for. 

The rates of change of fuel and emitter temperature and reactor power depend on 
the temperature coefficient of reactivity in the reactor. 

(1) Case 1. When fuel -emitter separation occurs in a single diode, the fuel tem- 
perature rises about 200° R (111 K) in the first 5 seconds. At this time, the emitter 
temperature reaches a minimum value about 600° R (330 K) below its initial value. The 
fuel and emitter temperature continue to increase, thereafter, until equilibrium is 
reached about 300 seconds after separation. The final fuel temperature is about 2500° R 
(1390 K) higher than its initial value while the final emitter temperature is the same as 
its initial value. H, as is indicated in reference 1, the maximum initial fuel tempera- 
ture is 4860° R (2700 K) and the UOg melting point is 5460° R (3030 K), the fuel will 
start to melt after the temperature has risen 600 u R (330 K). From figure 4, the fuel 
would start to melt about 17 seconds after the fuel separates from the emitter. When the 
fuel separates from the emitter in a single diode there would be no detectable change in 
reactor power because the reactor reactivity change associated with fuel-emitter separa- 
tion in a single diode is insignificant. 

(2) Case 2. When fuel-emitter separation occurs in all the fuel cells but the tem- 
perature coefficient of reactivity in the fuel is zero, reactivity and power will increase 
as the emitter temperature decreases because of the negative temperature coefficient of 
reactivity associated with the emitter. After about 5 seconds, the emitter temperature 
reaches a minimum and starts to increase and reactor reactivity decreases. About 

20 seconds after fuel-emitter separation the reactivity is too low to support a power rise, 
and reactor power decreases until it returns to the initial value about 300 seconds after 
separation. Both the fuel and emitter temperature reach maximum about 100 seconds 
after separation and decrease to equilibrium temperatures. The emitter equilibrium 
temperature is the same as the initial temperature and the fuel-equilibrium temperature 
is about 2500° R (1390 K) higher than the initial temperature. For the reference de’sign, 
the fuel would start to melt about 9. 5 seconds after separation. 

(3) Case 3. The effect of a positive temperature coefficient of reactivity (8. 0x10”^ 
$/°R) in the fuel was also investigated. In this case, the fuel temperature and power 
both rise faster than in the previous case because of the excess reactivity addition as- 
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sociated with the fuel temperature rise. For the reference reactor design conditions, 
the fuel would start to melt about 8 seconds after separation. 

Fuel-emitter separation in all fuel cells would be accompanied by reactor power 
increases, as indicated in figure 4. Separation could probably be detected by neutron 
flux monitors and reactor power could be reduced before fuel melting occurred in the 
reactor. In case 2, a 20 percent increase in power would occur within 1.5 seconds leav- 
ing about 8. 5 seconds to detect the abnormality and reduce power. Approximately 
7 seconds warning time would be available in case 3. 

It should be mentioned at this time that the 8. 5- and 7-second warning times are 
based on an arbitrarily selected 20 percent power rise and depend upon the operating 
temperature of the fuel, fuel- and material-temperature coefficients of reactivity, and 
on the number of fuel cells which experience separation. 

H the absolute values of the temperature coefficients of reactivity are smaller, or 
if only a few diodes experience fuel-emitter separation, the warning time would be re- 
duced because the rate of power rise would be reduced. It would, therefore, take longer 
to detect a given power rise. When only a few cells experience fuel-emitter separation, 
it would be difficult to prevent fuel melting without monitoring the temperature of all 
fuel cells. 

The calculated values of fuel temperature shown in figure 4 are conservative (high) 
because of several simplifying assumptions made in developing the mathematical model. 
Complete separation is assumed along the entire fuel-emitter interface. Actually some 
thermal contact would probably exist at some locations between the fuel and emitter after 
separation. Also, the thermal emissivities of the radiating surfaces were assumed to 
be equal (0. 2) and independent of time. The mass transfer between the fuel and emitter 
surfaces was also neglected. Actually, the emissivity of the uranium dioxide is temper- 
ature dependent (therefore, it is time dependent) and a substantial fraction (about 
20 percent is indicated in ref. 2) of the heat may be transferred across the fuel-emitter 
gap by mass transfer. This would reduce the thermal resistance between the fuel and 
emitter. The resulting change in fuel and emitter temperature and reactor power would 
be smaller. It would, therefore, take longer to detect a particular power change, but 
the fuel temperature change for a particular time after fuel-emitter separation would be 
less. The time -dependent effects of mass transfer have not been included in this study 
because the mass transfer cannot adequately be expressed as a linear function of tem- 
perature as required for solution by the AIROS (ref. 5) computer program. Estimates 
of the maximum fuel-temperature rise in a single fuel cell including both radiation and 
mass transfer indicate the fuel-temperature rise would still be greater than 1000° R 
(555 K) and the fuel would melt. 

The high fuel temperatures would not persist indefinitely. After fuel-emitter sepa- 
ration the hot fuel would vaporize and redeposit on the cooler emitter surface reducing 
the maximum fuel temperature. 
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METHODS OF REDUCING FUEL TEMPERATURE 


There are several potential methods of reducing the fuel temperature rise following 
fuel-emitter separation in an externally fueled diode. These include 

(1) Use a heat-conducting gas annulus between the fuel and emitter 

(2) Provide heat-conducting fins to transport the heat from the fuel to the emitter 

(3) Alter the material in the fuel region to produce a negative temperature coefficient 

of reactivity 

(4) Reduce the requirements on operating variables such as power density and cool- 

ant temperature 

The first two methods will reduce the fuel temperature whether the separation 
occurs in only a few or many diodes. The third method will reduce the fuel temperature 
only when sufficient diodes are affected such that reactor power will be reduced. The 
last method either results in a larger reactor, reduced electric power output or both and 
is, therefore, the least attractive of the four methods. Only the first two methods were 
investigated in detail. 


Gas Annulus 

The fuel-emitter separation problem might be avoided by intentionally introducing a 
gap containing a heat-conducting inert gas (such as helium) between the fuel and emitter. 

A slightly higher design fuel temperature would be accepted in order to avoid the large 
fuel-temperature rise resulting after fuel-emitter separation. Thermodynamically, a 
heat -conduction resistance (the gas space) would be substituted for a thermal-radiation 
resistance. Therefore, for a given heat flux, the conduction across the helium gap 
would result in a smaller temperature difference between the fuel and emitter than 
would be expected with radiant heat transfer between the fuel and emitter. 

The fuel-temperature rise resulting from variations in gap size would be about 
130° R per mil (28 000 K/cm) increase in gap thickness if helium is used in the gap. Dur- 
ing normal operation at 100 percent power in the reference system, an increase in gap 
size of about 5 mils (0. 013 cm) could be tolerated without melting the fuel. Without the 
gas annulus, even less than 1-mil (0. 003 -cm) separation between the fuel and emitter 
would be sufficient to cause the fuel to melt. The increase in normal operating fuel tem- 
perature caused by introducing the gas annulus might be acceptable. If not, either- the 
fuel thickness, power density, power distribution, or coolant temperature would have to 
be reduced to compensate for the excess temperature rise. 

One potential problem associated with this method, however, is that the fuel might 
be vaporized, transferred across the helium annulus, and redeposited on the emitter. If 
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sufficient fuel is redeposited on the emitter, the redeposited fuel might again separate 
from the emitter if some perturbation in the system caused the fuel temperature to rise. 
The rate of UOg vaporized from the fuel surface was approximated using the relation 
presented in reference 6 for the evaporation of a material into a vacuum, 


m = 5. 833X10 -2 P 



where 

m mass of UOg vaporized from fuel surface per unit time, g/(cm )(sec) 

M gram molecular weight of material 

P vapor pressure of material, mm Hg 

T temperature, K 

Using the procedure of Ward and Ruch (ref. 2), the net mass-transfer rate from the hot 
surface to the cooler surface is 


Am = 5.833X10" 2 




where the subscripts H and C refer to the hotter and cooler surfaces, respectively. 
The fuel vapor pressure P was approximated using the empirical correlation presented 
in reference 7 


log P = 13. 340 - 


3. 7337X10 4 3. 67X10 6 , 2.4638X10 9 

— . — - 4 - 4 


where P is in millimeters of mercury and temperature T is in K. 

The time required to transfer fuel across the helium gap is shown in figure 5 as a 
function of the relative thickness (r - Rg)/(R,j - Rg) of fuel transferred and helium gap 
width. Three modes of energy transfer were allowed in the calculation of the tempera- 
ture rise across the helium gap; conduction, radiation, and mass transfer. The calcu- 
lations are similar to those described in reference 2 except conduction is allowed across 
the helium. The results shown in figure 4 indicate that the fuel is transferred across the 
gap faster for large fuel -emitter gaps than for small gaps. This is expected because 
the larger gaps cause greater temperature differences across the gap and, thus, a 
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Figured. -Time required to transfer y percent of fuel 
thickness across helium gas gap. 


greater driving force for mass transfer. However, even with a small 0. 12 mil 
(3. 05x10 cm) helium gap, 90 percent of the fuel would be redeposited on the emitter 
after about 7000 hours of reactor operation and the fuel would again be subject to melt- 
ing caused by fuel-emitter separation if the helium does not migrate into the new fuel - 
emitter gap. 

The presence of helium in the fuel-emitter gap may impede the vapor transport of 
UOg across the gap. Gluyas (ref. 8) showed that, at a particular temperature, the pre- 
sence of argon at a pressure of 1 atmosphere could reduce the rate of mass transfer of 
UOg by an order of magnitude. However, since the rate of mass transfer of UOg across 
a helium gap is questionable, and because of the difficulty associated with maintaining 
any prescribed gap between the fuel and emitter during reactor operation, this method 
of avoiding fuel melting would probably have to be demonstrated. 


15 



Heat Conducting Fins 


Another potential method of avoiding fuel melting is to provide fins which would 
conduct heat from the fuel to the emitter and, thus, reduce the maximum fuel tempera- 
ture after the fuel separates from the emitter. A sketch of one possible finned therm- 
ionic diode concept is shown in figure 6. Radial-tungsten fins are incorporated within 
the fuel. The thermal conductivity of the tungsten is about 60 times better than that of 
the fuel. If fuel-emitter separation occurs, heat can be conducted from the fuel into the 

' V P ' ' s 

fin, and radially conducted along the fin into the tungsten emitter. Any thermal expan- 
sion induced in the fuel will tend to improve the contact between the fuel and the fin. 

The fuel and reactor diameters must be increased in order to add sufficient fuel to com - 
pensate for the additional fin material and to maintain nuclear criticality. The maxi- 
mum temperature in the fuel was conservatively estimated by assuming only azimuthal 
heat conduction in the fuel, such that all of the heat generated in the fuel is radially 
transferred along the tungsten fin to the emitter. Reactor power and the emitter heat 
flux were assumed to be 208 kilowatts electric and 68 Btu per second per square foot, re- 
spectively (see pp. 40 and 41 of ref. 1). The reactor length was maintained at 18 inches 
(46 cm). Neutron and gamma heating in the tungsten fin is insignificant compared with the 
total heat generated in the fuel and was neglected. To simplify the mathematical analysis, 
the fin width was assumed to be proportional to the radius as shown in the sketch of a typi- 
cal section of the fuel region shown in figure 7. The planes defined by 6 = 0 and 6 = ©g 
are surfaces of symmetry, and the outer fuel surface is an adiabatic surface. Based 
upon these boundary conditions and assuming that the temperature of the fin root is equal 


Coplant 
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Fuel 


/ 



to the emitter temperature, the maximum fuel temperature can be approximated by the 
following equation which is derived in appendix B: 



= T 

max E 



''t 



j 


With constant fuel and emitter thermal conductivities (kp and kg), the maximum 
fuel temperature depends on the volume fraction of fin in the fuel region 9 9 2 > the 

number of fins n, the fuel power density q F , and the outer radius of the fuel R^. 

The power density and fuel outer radius required to produce an emitter heat flux of 
68 Btu per second per square foot (77 W/cm ) and to maintain nuclear criticality de- 
pend upon the volume fraction of fin in the fuel. Therefore, q F and Rg were calculated 
for various volume fractions of fin in the fuel. The temperature difference (Tp max - 
Tp) is shown in figure 8 as a function of the volume fraction of fin in the fuel and the 
number of fins. As more fins are added, the maximum fuel temperature is reduced be- 
cause the azimuthal heat conduction path in the fuel is reduced. As the volume fraction 
of fins in the fuel region is increased (each fin is made thicker), the fuel temperature is 
initially reduced because the heat transfer area of the fin is increased. However,, the 
height of the fin must also be increased as the width is increased in order to add suffi- 
cient fuel to maintain reactor criticality. When the increase in fin height exceeds the 
increase in width, fuel temperature will have reached a minimum and will increase 
thereafter. The minimum fuel temperatures are obtained with fin volume fractions be- 
tween 0. 80 and 1. 0. 
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0 .2 A .6 .8 1.0 

Fraction of fuel volume displaced by fin, 

Figure 8, - Effect of fin design on maximum 


fuel temperature. Emitter heat flux, 77 watts 
per square centimeter (68 Btu/(sec)(ft 2 )). 

These results are shown in a different form in figure 9. The number of fins required 
to obtain various temperature rises is shown as a function of the volume fraction of fin 
in the fuel region. The relative volume of the reactor and the relative outer radius of 
the fuel required for criticality is also shown as a function of the volume fraction of fin 
in the fuel region. Although the lowest fuel temperatures would be obtained with fin 
volume fractions greater than 80 percent (fig. 8), the large increase in reactor size re- 
quired for criticality would be impractical for space power considerations. In the reac- 
tor design of reference 1 the value of T-r, - T„ corresponding to fuel melting is 

q -C j iXl3*X Jui 

I860 R (1030 K). If fuel melting is to be avoided, the appropriate number of fins must 
be selected such that the temperature difference is less than 1860° R (1030 K), The 
number of fins must be selected such that a compromise is reached between final reactor 
volume and ease of fabricability of the fin. For example, four fins require a fin volume 
fraction of 0. 48 in order to avoid fuel melting as shown in figure 8. However, 

as indicated in figure 9, the reactor volume would have to be increased by a factor of 
2. 5 to maintain criticality. The selection of 10 fins requires a fin volume fraction of 
0. 02. Although this requires a reactor volume increase of only 2 percent (from fig. 7), 
the fin root width w required with 10 fins (n = 10) and a volume fraction of 0. 02 is 
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Figure 9. - Number of fins required and ef- 
fect of fin size on reactor volume. 


2vR a — 
o o 


W = 


n 


= — ^°- 266 )f°- ° 2 ) = 0. 0033 inch (0. 0084 cm) 

10 

It might be difficult to fabricate and assemble a fueled diode having a tungsten fin 
3. 3 mils (0. 0084 cm) wide, 0. 5 inch (1. 27 cm) high and about 1 inch (2. 54 cm) long. 

For the reference design externally fueled thermionic diode, eight fins having a fin 
volume fraction of 0. 03 and a fin root width of 0. 006 inch (0. 0152 cm) would keep the 
maximum fuel temperature below the melting point. The reactor volume and fuel outer 
radius would have to be increased by only about 3 percent and 1. 7 percent, respectively, 
to maintain criticality. 


CONCLUSIONS 


The temperature response following separation of the fuel from the emitter in an 
externally fueled thermionic diode was calculated neglecting the mass transfer across 
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the fuel-emitter gap. The results of the study on the reference design externally fueled 
thermionic diode indicate 

1. If fuel -emitter separation occurs in many diodes, the fuel would start to melt 
within about 8 to 10 seconds after separation if reactor power is not reduced. The power 
rise associated with this accident is sufficiently large to be detected by neutron flux 
monitors and could allow reactor power to be reduced before the fuel melts. 

2. If, however, fuel-emitter separation occurs in only a few diodes, the fuel in 
these diodes could melt about 17 seconds after separation and there would be no detect- 
able increase in neutron flux. Unless the temperature of the fuel material in each diode 
is monitored, there may be no way of detecting separation in time to prevent fuel melt- 
ing. 

Neglecting the effect of mass transfer across the fuel-emitter gap results in con- 
servative (too high) fuel temperatures. However, steady-state calculations, including 
both the effects of mass transfer and radiation, indicate that the maximum fuel temperature 
would be less but fuel melting would still occur. The reduction in thermal resistance 
caused by the inclusion of mass transfer would also result in slower response times for 
fuel and emitter temperature and reactor power. 

Several methods of reducing the fuel temperature rise after separation were inves- 
tigated. The most attractive method required the redesign of the diode to include radial 
tungsten cooling fins in the fuel region which would allow conduction of heat from the 
fuel through the fin and into the emitter. It is desirable to minimize total fin volume 
because, to maintain nuclear criticality, reactor volume must be increased to compensate 
for the addition of fins. Therefore, a compromise must be reached between reactor 
volume and fabrication capability. For the reference design thermionic diode, fuel 
melting could be avoided by inserting eight fins in a diode. The required fin width is 
about 6 mils (0. 015 cm). The reactor volume and fuel outer radius would have to be 
increased about 3 percent and 1. 7 percent, respectively. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 8, 1969, 

120-27-05-20-22. 
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APPENDIX A 


MATHEMATICAL MODEL 

The time -dependent temperature distribution in the externally fueled thermionic 
diode was calculated using the digital computer program AIROS (ref. 5). This program 
allows the solution of simultaneous linear differential equations of first order and first 
degree. These equations are coupled to the space -independent reactor neutronics equa- 
tion. The heat transfer through the diode was approximated using a one -dimensional 
(radial), time -dependent model. The model allows heat conduction through all compo- 
nents of the diode. Heat is removed by forced convection to an incompressible liquid 
coolant flowing through the coolant tube. Radiation between the fuel and clad after sepa- 
ration is approximated by a linear function of the temperature difference between the 
fuel and clad. The diode efficiency r\ is assumed constant throughout the transient. 

The following equations represent a heat balance over each major diode component: 

Coolant: 

dT R 

P B C P b Vb = U P A P (T P " T B> “ 2F B C p B < T B “ T in> < A1 > 

Coolant tube: 

dT 

Pp C Pp V p - Wi - V - Vp (t p - T B> <a2) 

Insulator: 

dT, 

Pi °pi Vi ~i = UcA c (Tc ' Ti> ' UiAi<Ti ' v (A3) 

Collector: 

dT p 

pc c p c v c IT = (1 ' ^E^E " T C> “ u c A c< T c - T i) < A4 > 

Emitter: 

dT,-, 

Pe c p e v e --f = W T F,s - t e) - Wi - V + Wbe^bp - t e) < A5 > 
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Fuel: 


dT„ 

P F C P f V F _ ^' = q F V F " U F A F^ T F,s " T E> “ U BF A BF^ T F " T BP^ ( A6 ^ 

Unless otherwise indicated, the overall conductances t^A., which are constant 
were calculated for initial, steady -state conditions. The bulk coolant temperature T B 
is 



T. + T 
m o 


(A7) 


All component temperatures except the fuel surface temperature T p g are taken at the 
radial midpoint of the component. The overall conductances used to describe the heat 
transfer between the pipe and coolant U A , insulator and pipe U T A T , and collector 
and insulator U^A^, are 


U A = 
p p 


2 TTl 


In I 


+ R2^ 

2Ri i 


R t h 


U I A I = 


2 TTl 


In 


2R r 


\ R 1 + R 2 , 


In 


( 2R 0 j 


U C A C ~ 


2 Til 


In 


^ R 3 + R 4^ 

, 2Ro I 


In 


/ 2R 3 N 
V R 2 + R 3/ 


(A8) 


(A9) 


(A10) 


The overall conductance between the emitter and collector includes several modes 
of heat transfer. This conductance was obtained from a typical diode performance curve 
(ref. 1) by empirically expressing the heat flux as a function of temperature drop across 
the interelectrode gap for constant emitter current density. This relation is presented 
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in figure 10. If the accident occurred in a single cell, the current density would nearly 
be constant because all cells in a thermionic fuel assembly are connected in series. The 
empirical correlation for the overall conductance is 

U e A e = 0.20229 - 1.2186X10" 4 (T e - T c ) (All) 

where U E A E is in units of watts per °K, and T £ and are in units of °K. Using 
initial values of T £ (2000 K) and T c (1400 K), 


U e A e = 0. 244 


(hr)(°R) 


0. 129 


The overall conductance between the fuel and emitter after separation includes only 
radiant heat transfer and is expressed as 




2l7<7i(R 6 + AR) (t| iS + t| jS T e + T Fi5 t| + t|) 


R. + AR 



Figure 10. - Effect of input power on 
emitter temperature. Collector 
temperature, 2520° R (1400 K). 
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The final steady -state values of Tp g and Tp obtained using the initial fuel power den- 
sity qp were used to calculate UpAp. Prior to separation, the fuel and emitter are in 
thermal contact and Tp g « Tp. The thermal driving force in equations (A5) and (A 6) 
is t f“ t e instead of Tp g - Tp and 


UpAp - 


2 7 Tl 


In 



+ 



(A13) 


The difference between the fuel temperature Tp and the fuel-surface temperature 
Tp is assumed to be independent of time. 

Heat conduction from the fuel, through the end clad, and into the emitter is accounted 
for by the terms u bE A BE^ T BP “ T E^ and ^BF^BF^F ” T B p) in equations (A5) and (A6). 


U BE A BE 



(A 14) 


U BF A BF “ 


k F( R ?- R e) 

The temperature response of the end clad is calculated by 
dT 

P BP C p Bp V BP “ -jp = U BF A BF^ T F ” T BP^ " U BE A BE^ T BP ‘ T E> ( A16 > 

** 

* The reactor power and fuel assembly power density qp were approximated by the 
monoenergetic space -independent neutronics equation including six delayed neutron 
groups: 


2lt 


In 


rRg + R/y + 2AR\ 
y 2Rg + AR J 


(A 15) 
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(A 17) 


dqj, 

dt 


§_ 

L 


|[r(t) - 1] q p 



E 



dq ni ft* 

__I=J:q -Aq (A 18) 

dt L * 1 m 

6 

(3= 2^ h (A19> 

i = 1 

The effect of changes in the component temperatures on reactivity T(t) is approxi- 
mated using temperature coefficients of reactivity for each of the j components as 
follows 


r(t) = r(0) C.[T.(t) - T.(0)] 

j 


(A20) 


where the reactivity r is in dollars. 

Equations (Al) to (A 19) are solved on a digital Computer to determine the component 
temperature and power response as a function of time. 

Calculations were performed using the cell dimensions presented in table I, fuel 
element operating conditions given in table II, initial temperatures shown in table HI, 
and component property values shown in table IV. 
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TABLE IV, - COMPONENT PROPERTIES 


Component 

Material 

Property 

Value 

Refer- 





ence 

Coolant 

Li 

p, lb/ft 3 ; g/cm 3 

26.2; 0.420 

9 



Cp, Btu/(lb)(°R); J/(kg)(K) 
k, Btu/(sec)(ft)(°R); W/ (cm) (K) 

0.996; 4169 

9 



1.04X10' 2 ; 0.648 

9 

Coolant tube 

Nb 

P, lb/ft 3 ; g/cm 3 

535; 8.57 

10 



C p , Btu/(lb)(°R); J/(kg)(K) 

0.085; 355 

10 



k, Btu/(sec)(ft)(°R); W/(cm)(K) 

9. 7xl0~ 3 ; 0. 61 

10 

Insulator 

A1 2°3 

p, lb/ft 3 ; g/cm 3 

240; 3. 85 

11 



C p , Btu/(lb)(°R); J/(kg)(K) 

0.306; 1280 

11 



k, Btu/(sec)(ft)(°R); W/(cm)(K) 

9.7X10' 4 ; 0.056 

11 

Collector 

Nb 

p, lb/ft 3 ; g/cm 3 

535; 8.57 

10 



C p , Btu/(lb)(°R); J/(kg)(K) 

0.085; 355 

10 



k, Btu/(sec)(ft)(°H); W/(cm)(K) 

9. 7X1 0'"' 3 ; 0.61 

10 

Emitter 

W 

p, lb/ft 3 ; g/cm 3 

1210; 19.4 

4 



C p , Btu/(lb)(°R); J/(kg)(K) 

0.035; 150 

4 



k, Btu/(sec)(ft)(°R); W/(cm)(K) 

1.78X10" 2 ; 1. 11 

4 

Fuel 

uo 2 

p, lb/ft 3 ; g/cm 3 

684; 11.0 

3 


1 

C p , Btu/(Ib)(°R); J/(kg)(K) 

0.08; 300 

3 



k, Btu/(sec)(ft)(°R); W/(cm)(K) 

2. 78x10 " 4 ; 0.0174 

3 

Clad 

w 

p, lb/ft 3 ; g/cm 3 

1210; 19.4 

4 



C p , Btu/(lb)(°R); J/(kg)(K) 

0.035; 150 

4 



k, Btu/(sec)(ft)(°R); W/(cm)(K)j 

1.78X10" 2 ; 1.11 

4 
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APPENDIX B 


ESTIMATED HEAT TRANSFER FROM EXTERNALLY FUELED 
THERMIONIC DIODE CONTAINING RADIAL FINS 

The temperature distribution in an externally fueled diode containing radial-cooling 
fins connecting the fuel with the emitter was estimated assuming only azimuthal heat 
conduction in the fuel and only radial heat conduction in the fin. Using the sketch shown 
in figure 7, the temperature distribution in the fuel region can be conservatively esti- 
mated if radial and axial heat transfer are neglected after fuel -emitter separation occurs. 
The Fourier equation is used to approximate the azimuthal heat conduction in the fuel 
having a fission (power) density of q^: 

J_ d^T _ _ ^F 
r 2 d0 2 k F 

If is taken as the plane of symmetry between adjacent fuel regions and the tempera- 
ture at the fuel -fin interface is T B> the boundary conditions are 


— = 0 at 0 = 0<j 
d0 1 

T = T b at 0 - 0j 

The azimuthal temperature distribution is: 

t=t b + ^- [e(20 2 - e) - 0 i(20 2 - &1 )] 

The steady -state heat transfer in the fin can be approximated with a heat balance 
over an incremental radius of fin, dr; that is, the heat transferred into the fin per unit 
length of the fuel is equal to the heat conducted radially out of the fin per unit length of 
fin: 


r< ^F^2 


- 0 j) dr - “k E — 
dr 



dr 
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Reduced to the simplest algebraic form, this equation becomes 

d_ / dT\ 
dr \ dr/ k E 

If no heat transfer is assumed from the outer radius of the fin and the tempera- 
ture at the fin root Rg is assumed to be the same as the emitter temperature, the 
boundary conditions for the fin are 



— =0 at r = R„ 
dr 7 

T = T e at r = Rg 

The solution to the fin heat transfer equation is 



The maximum temperature difference between fuel and emitter is 



’ The number of fins in a diode n is equal to the circumference of the fuel 2jrr 
times the fin volume fraction divided by the width of each fin 20^r: 
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2irr — 


20 l r 


or 


n = 


IT 


e 


2 


The maximum temperature rise across fuel and fins can be expressed in terms of the 
fuel volumetric heating rate q^, material thermal conductivity ratio kp/k^, fuel 
dimensions Rg and R^, fin volume fraction 0 j/flgj an< ^ number of fins connecting the 
fuel with emitter n: 
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